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Background Neutrino oscillation probabilities, which are being measured in long-baseline experiments, depend on neutrino 

energy. The energy in a neutrino beam, however, is broadly smeared so that the neutrino energy in a particular event is 

not directly known, but must be reconstructed from final state properties. 
Purpose We investigate the contributions from different reaction mechanisms on the energy-reconstruction method widely 

used in long-baseline neutrino experiments. We clarify how the oscillation signal depends on the energy reconstruction 

method. 

Methods The Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) transport model is used for a detailed study of neutrino-nucleus 
events. 

Results Difference between the true-QE and QE-like cross sections in MiniBooNE experiment is investigated in details. It 
is shown, that fake QE-like events lead to significant distortions in neutrino energy reconstruction. Flux-folded and 
unfolded cross sections for QE-like scattering are calculated as functions of both true and reconstructed energies. Flux- 
folded momentum transfer distributions are calculated as functions of both true and reconstructed momentum transfer. 
Distributions versus reconstructed values are compared with the experimental data. Also presented are the conditional 
probability densities of finding a true energy for a given reconstructed energy. It is shown, how the energy reconstruction 
procedure influences the measurement of oscillation parameters in T2K experiment. 

Conclusions For the reconstruction procedure based on quasielastic (QE) kinematics, all other reaction channels beside true- 
QE scattering show a shift of the reconstructed energy towards lower values as compared to the true energy. On average 
in the MiniBooNE and T2K experiments the shift is 100 - 200 MeV and depends on energy. The oscillation signals are 
similarly affected. These uncertainties may limit the extraction of a CP violating phase from an oscillation result. 

PACS numbers: 



I. INTRODUCTION 

The extraction of neutrino oscillation parameters nec- 
essarily requires the knowledge of the neutrino energy. 
Since neutrino beams are always produced as secondary 
decay products their energy is not sharp, but widely dis- 
tributed. The energy of the incoming neutrino in a given 
event thus has to be reconstructed from observed proper- 
ties of the final state. Often QE scattering is used for this 
reconstruction, because for this process the incoming en- 
ergy can be uniquely inferred from an observation of just 
the outgoing lepton if the target is a nucleon at rest. This 
becomes considerably more complicated, however, when 



nuclear targets are used [TJ [5] . For these, the difficulty 
is twofold. First, the reaction process of QE scatter- 
ing must be unequivocally identified. Second, even then 
nuclear effects can smear out the reconstructed energy. 
The first difficulty is the more serious one, since other 
reaction mechanisms may look indistinguishable for the 
experiment. Furthermore, final state interactions make 
it very difficult to identify the initial QE scattering on 
a bound, Fermi-moving nucleon inside the nucleus. This 
is even more so if no outgoing nucleons are observed, as 
is at present the case in all experiments using Cerenkov 
detectors. 

Recently quasielastic (QE) neutrino scattering has at- 
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tracted a lot of attention because the cross sections 
reported lately [21 0] by the MiniBooNE collaboration 
at E v < 2 GeV are 30 - 40% higher than the cross 
section measured by old bubble-chamber experiments 
( ANL ,BNL ,FNAL , CERN,IHEP) in the 70s-80s. These 
higher data could be well described only by assuming 
an axial mass that was significantly higher than the 
world average value of Ma ~ 1.0 GeV [SH3- The sit- 
uation is complicated by another modern experiment - 
NOMAD [5] - which at E v > 5 GeV reports M A = 
1.05 ± 0.08 GeV and cross sections which are in agree- 
ment with the old measurements. The Minerva experi- 
ment, operating with a higher energy neutrino flux, also 
does not derive high a Ma from their results [9]. 

At first sight, the modern experiments have the advan- 
tage of huge statistics with millions of events recorded. 
The complication arise, however, from the fact that these 
experiments all use nuclei as targets. All the measure- 
ments are thus influenced by nuclear effects. These affect 
primarily the event identification and with it the energy 
reconstruction. This is crucial for the extraction of the 
cross section as a function of neutrino energy [7] and in- 
fluences the extraction of oscillation parameters |10j . 

The true charged current (CC) QE scattering is defined 
as neutrino scattering on bound neutron in the nucleus, 
resulting in a muon and a proton 

true QE: vn — > jj,~p 

Such events can be identified quite well in a tracking 
detector, but are impossible to identify in a Cherenkov 
detector, because it is "blind" to outgoing neutrons and 
low-energy outgoing protons. Thus, in a Cherenkov de- 
tector (e.g. MiniBooNE and T2K), the signal is defined as 
a single Cherenkov ring from the outgoing muon, which 
can also be identified by its decay electron. No further 
rings, possibly originating from pions or other mesons, 
should appear in such an event. Thus, CC QE-likc events 
are defined as those with 1 muon, mesons and any num- 
ber of nucleons in the final state. 



Not all events that look like QE in the Cherenkov de- 
tector have a true QE origin. In the MiniBooNE and T2K 
energy regimes one complication arises from processes in 
which a pion or a A resonance are produced in the ini- 
tial neutrino vertex, e.g., up — > [i~pir + or vp — > A ++ , 
and then are absorbed in the nucleus through final state 
interactions. Thus the pion is not seen in the final state 
(so-called 'stuck-pion' event) and such event is counted 
as QE-like even though it is of not true-QE origin. In an 
earlier publication [TT] we have shown that there is indeed 
a strong entanglement of QE scattering and pion produc- 
tion events and that, consequently, Cherenkov detectors 
always see too high a cross section for QE. A second com- 
plication is the presence of multi-nucleon events in which 
the incoming neutrino interacts with e.g. 2 nucleons (so- 
called 2p-2h events). Both of these complications lead 
to the so-called fake QE events. The "measured" QE- 
like cross section is then contaminated by these fake QE 
events. On the other hand, a nucleon produced in the ini- 
tial QE vertex may rescatter in the nucleus and produce 
a pion in the final state. This event will be disregarded as 
QE-like, even though it originates in a true-QE scattering 
process. 

The experimental groups try to account for nuclear ef- 
fects by using event-generators. Thus, the final 'data' 
for QE scattering invariably contain some model depen- 
dence and may suffer from imperfections in the models 
used. For example, the widely discussed large values for 
the axial mass (M A w 1.3 GeV) obtained by the Mini- 
BooNE collaboration are nowadays believed to be due to 
2p-2h excitations (see [HI [T3] and refs. therein). These 
excitations were not contained in the generator used and 
could thus not be removed from the data set. In addition, 
one has to be aware that in the MiniBooNE data shown 
in Fig. [T] neither the QE cross section itself nor the neu- 
trino energy have been directly measured. Instead only 
the flux-averaged double-differential (with respect to out- 
going lepton variables) QE-like cross section can directly 
be measured. Even the energy-separated QE-like cross 
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section (full squares in Fig. [I] denoted as "measured") 
as a function of neutrino energy is model-dependent; the 
same is naturally true for the QE-extracted cross sec- 
tion (open triangles in Fig. [I] denoted as "extracted"). 
In both cases the energy-dependence has been recon- 
structed. The stuck-pion processes should account for 
most of the difference between these two data sets. 
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FIG. 1: (Color online) The green squared points give the 
"measured" QE-like cross section obtained by the MiniBooNE 
collaboration [TJ. Subtracting from this the "QE-like back- 
ground" (see text), the "extracted" cross section shown by 
the red triangular points is obtained. All data are plotted vs 
reconstructed energy. The solid (red) line gives the result for 
the true-QE cross section, obtained in a calculation using the 
world average value for the axial mass of Ma = 1 GeV. 

It is the purpose of this paper to explore the full QE- 
like cross section that includes the events from various 
origins. These are clarified in sec. [TTJ. Special attention 
is given to the 2p-2h processes, which are discussed in 
sec. II C In sec. Ill we consider the sensitivity of the en- 



ergy reconstruction method used by Cerenkov detectors 
such as MiniBooNE to the events of each origin in a more 
comprehensive way than done so far. In particular, we 
explore the dependence of the energy reconstruction not 
just on 2p-2h and A events, but also on QE-like con- 
tributions from pion background, higher resonances and 
DIS. Reconstruction of momentum transfer is discussed 



in sec. IV The influence of energy reconstruction proce- 



dure on measurement of oscillation parameters in T2K 
experiment is illustrated in sec. [V] At the end we sum- 
marize our findings. 



II. QE-LIKE CROSS SECTIONS 
A. Theoretical Method 

For our investigations of nuclear effects on the QE- 
like and extracted QE cross sections we use the Giessen 
Boltzmann-Uehling-Uhlenbeck transport model GiBUU. 
Its theoretical foundation and details of the practical im- 
plementation are described in detail in Ref. [14] . In this 
model we solve the approximate Kadanoff-Baym equa- 
tions for the time-development of the spectral phase 
space distributions of nucleons and mesons after a nu- 
cleus has interacted with an incoming particle. The 
model, which has been widely tested on a variety of nu- 
clear reactions, is an event generator that delivers at 
the end, i.e. after all final state interactions, the four- 
momenta of all asymptotically free on-shell particles. It 
contains all the relevant processes: QE scattering, back- 
ground pion production and nucleon resonance dynam- 
ics are included; all resonance properties are taken ei- 
ther from the PDG listings or (for electromagnetic vec- 
tor couplings) from the MAID analysis [T51 [TB]. The 
model has recently been extended by also including ini- 
tial two-particle- two-hole (2p-2h) interactions [17]. At 
higher energies pQCD expressions encoded in the high- 
energy event generator PYTHIA are used [H] . The model 
thus contains all the relevant processes for the initial neu- 
trino interaction. All calculations reported in this paper 
have been performed with an axial mass Ma — 1 GeV 
for QE scattering; model II from Ref. [17 is used for the 
2p-2h contributions. Once the particles are produced in 
the initial vertex they are being transported out of the 
nucleus. During this process elastic and inelastic scat- 
tering can take place, including charge transfer, particle 
production and absorption. 



4 



B. True QE and QE-like cross sections 
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FIG. 2: (Color online) True QE and QE-like cross sections for 
neutrino scattering off carbon. Model II from [17] has been 
used for the 2p-2h interactions. 

Fig. [2] shows true-QE (solid curve) and QE-like (short- 
dashed curve) cross sections versus neutrino energy as 
obtained in the GiBUU model for neutrino CC scatter- 
ing off carbon. The QE-like events have been identified 
using the same method as that employed with Cerenkov 
detectors. As one can see from Fig.[2j the Cherenkov de- 
tector sees almost all true QE events (dash-dotted curve) , 
but also a large part of the fake QE-events. Such events 
are called "QE-like background" by MiniBooNE. So the 
resulting QE-like cross section is about 30-40% higher 
than the true-QE one. The fake events, as shown in 
Fig. [2] may come from A and other resonances decay- 
ing to mesons, 1-pion background and, at higher ener- 
gies, from DIS when the outgoing mesons are absorbed 
in the nucleus during the FSI. In addition, there may be 
events present in which the initial interaction involved 
more than one nucleon (the so-called 2p-2h processes, 



see discussion in sec. II C ) 



In MiniBooNE, the QE-like background is not directly 
measured, but obtained from the NUANCE neutrino 
event generator which in turn was adjusted to the mea- 
sured pion yields [7J; this generator contained no 2p-2h 
processes. This background should correspond to the 



sum of fake events originating from resonances, 1-pion 
background, and DIS in the GiBUU. Subtracting QE-like 
background from the "measured" QE-like cross section, 
MiniBooNE presented its QE result (which is meant to be 
sensitive to the axial mass) as "extracted" cross section 
(open triangles in Fig. [IJ . 

The neutrino energy for each event has been recon- 
structed. The latter is necessary because in any exper- 
iment the neutrino beam involves broad energy distri- 
bution, and the true energy for a given event is thus 
not known. In reporting the inclusive cross sections (es- 
pecially at high neutrino energy), the experiments can 
mostly rely on calorimetric measurements of the energy 
of final hadrons (MINOS, NOMAD). In contrast, for a 
specific reaction mechanisms, such as QE scattering, the 
energy has to be reconstructed from muon observables 
only. The latter can lead to inaccuracies in the recon- 
structed energies that are larger than previously assumed 
if the actual reaction process is not correctly identified. 



C. Influence of 2p-2h interactions 

As we see from Fig. [TJ removing all the "stuck-pion" 
QE-like events does not lead to an agreement with the 
conventional theoretical calculation of true-QE cross sec- 
tion. The "extracted" cross section still is considerably 
higher than calculations using the world average axial 
mass of 1 GeV. The remaining difference can be at- 
tributed to a significant amount of non-QE many-body 
excitations in the QE-like cross section; for a compact 
review of other possible explanation see the introduction 
in [17] • Since the NUANCE generator used by Mini- 
BooNE does not contain 2p-2h excitations, they are not 
subtracted as a part of the QE-like background, and thus 
are still present in the "measured" cross section. 

By combining the random phase approximation (RPA) 
with a calculation of 2p-2h contributions Martini et al. 
have obtained a good description of the MiniBooNE data 
QH HU HD]. As expected, the RPA correlations have 
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most effect at forward angles where the squared four- 
momentum transfer Q 2 to the nucleus is small. They 
die out with increasing angle and with decreasing muon 
energy, i.e. increasing energy transfer. These results have 
been confirmed by detailed calculations of the Valencia 
group [13l [21] and the effects of such interactions have 
been explored in detail in Ref. [17] (see also Ref . [22] ) ■ 

Since all measurements involve a flux-average over a 
usually broad band of incoming neutrino energies we 
have developed a short-cut to the full theoretical treat- 
ment of such 2p-2h processes by parametrizing a flux- 
averaged matrix element. This method has allowed us 
to implement these processes in an actual event genera- 
tor, GiBUU, and to investigate observable consequences 
of these initial 2p-2h interactions mainly on knock-out 
nucleons. In [17j we have shown that the measured 
double-differential cross sections (corrected for stuck- 
pion events) could be described quite well so that our 
method seems to contain the relevant physics. There we 
had used a model for the 2p-2h part of the hadron tensor 
that consisted of the transverse projector, thus modeling 
the transverse character of 2p-2h interactions, modified 
by an explicit energy-dependence (model II in [17)). 

The 2p-2h cross section leading to fake QE-like events 
is also shown in Fig. [2] It is seen that at low energies, 
E v < 0.4 GeV, fake events come only from 2p-2h con- 
tributions because at these low energies processes involv- 
ing resonance excitations or pions are kinematically sup- 
pressed. Above this energy the main contribution comes 
from the 2p-2h mechanism and A excitation followed by 
its decay. With increasing energy the contribution of 
fake events grows, because new channels (production of 
higher baryonic resonances, DIS) open up. The relative 
contributions of fake events of 2p-2h (dash-dotted line) 
and stuck-pion (dashed line) origins are shown in Fig. [3] 
While the 2p-2h contribution is quite flat at about 20- 
30%, the pion contribution rises with energy, reflecting 
the threshold for pion production. 

In |17j we had determined the matrix element such that 
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FIG. 3: (Color online) Ratios of fake QE-like cross sections 
of 2p-2h and 'stuck pion' origins to the true-QE one. 



10 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

E v (GeV) 

FIG. 4: (Color online) QE-like cross section originating from 
QE and 2p-2h processes only (solid line) and from all pro- 
cesses (dashed line) within the GiBUU calculations. "Mea- 
sured" and "extracted" MiniBooNE data points are the same 
as in Fig. [I] The difference between them (open circles) is 
to be compared with the GiBUU 'stuck pion' cross section 
(dotted line). All data are plotted vs reconstructed energy, 
whereas the theoretical curves are plotted vs true neutrino 
energy. 

the sum of true-QE and 2p-2h contributions fitted the 
"extracted" MiniBooNE data. This is shown in Fig. [4] 
where the solid ("true-QE + 2p2h") line is the GiBUU 
model calculation that includes only true-QE and 2p-2h 
cross sections. Even with this fit, as illustrated in Fig. [4] 
the "measured" data points still do not agree with our 



curve for the total QE-like cross section. The latter is 
shown by the dashed ("all") line and includes all pro- 
cesses that lead to QE-like final state. 

As shown in the previous section, the absolute contri- 
bution of fake 'stuck-pion' QE-like events (that is, the dif- 
ference between the dashed and the solid curves in Fig. |1J 
also shown as dotted curve) is zero for E v < 0.4 GeV and 
slowly grows with increasing energy. The MiniBooNE re- 
sults (open circles), however, show quite a different pic- 
ture. The contribution of fake events is largest at low 
energies and decreases further as energy grows (open cir- 
cles in Fig. [I]). The theoretical "all" and "true-QE + 
2p2h" curves do not agree with the data, both have a 
noticeably different shape. As we will show later in this 
paper, the resolution of this seeming contradiction lies 
in the fact that in Fig. [4] the data are plotted versus re- 
constructed energy whereas the calculated curves are all 
plotted vs true energy 



III. ENERGY RECONSTRUCTION BASED ON 
QE KINEMATICS 

To resolve contradiction shown in Fig.|4j let us consider 
the energy reconstruction procedure used by MiniBooNE 
and its influence not just on the QE, but also on the QE- 
like cross sections. As was shown already in [TTl |2"31 123] , 
a 2p-2h interaction, when leading to a final state with 
pions and thus recorded as QE-like event, is on average 
recorded with a reconstructed energy lower than the true 
energy. 

For QE scattering on a nucleon at rest the incoming 
neutrino energy is directly linked to the kinematics of 
the outgoing lepton and is thus known when lepton an- 
gle 9^ and energy are measured. Therefore, the for- 
mula used by MiniBooNE for the energy reconstruction 
is based on the assumption of QE scattering on a nucleon 
at rest [35] even though nuclear targets with binding and 
Fermi-motion are used. The reconstructed (rec) neutrino 



energy is defined as 



ET = 



2(M„ - E B )E ll - {E 2 B - 2M n E B 



AM 2 



M n — Eb — En 



\kf_tl cost 



Here M„ the mass of the neutron, AM 



M, 



(1) 



and |fcp,| = ^JE 2 —m 2 l is the absolute value of the 3- 
momentum of the outgoing muon. This formula, there- 
fore, neglects any Fermi-motion effects; binding is taken 
into account only by a constant removal energy Eb > 0. 
It is essential to realize that use of this formula is justified 
only if the reaction mechanism has been identified as be- 
ing true QE scattering; admixture of any other reaction 
modes leads to an incorrect reconstruction of energy. In 
the following we will explore how large these errors actu- 
ally are. 

In a numerical simulation, e.g. with the MiniBooNE 
flux, the distribution of true neutrino energies is known. 
In GiBUU a flux-weighted Monte-Carlo sampling of true 
energies is then used to generate for each of them one 
event; from this event a reconstructed energy is calcu- 
lated according to Eq. This method thus directly 
corresponds to the method used in the experiment. Us- 
ing a sufficiently large number of true energies eventually 
the whole (E tIUC , E rcc ) plane is filled. The distribution of 
points in this plane is denoted by N (E tIUC , E ICC ) , where 
N is a 2-D density of the flux-averaged cross section. 1 
Integrating N over E truc and E rcc gives 

J N(E tIue ,E rec )dE Iec dE time 

= J<l>(E tro *)iTo«(E tim )dE trm = M (2) 

which is just the flux averaged cross section. 

The two-dimensional density N (E tTUC , E Tec ) for the 
MiniBooNE flux is shown in Fig. [5j For true-QE 
events the distribution is nearly symmetric about the 
line El ruc — El cc , the broadening comes from the Fermi 



1 In [24] the same value is denoted as d ^ (E; E Tec )$(E), see 
Eq.(7) there. 
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FIG. 5: (Color online) 2-D density of the MiniBooNE flux- 
averaged QE-like cross section N (E tTUC , E TCC ) versus true and 
reconstructed neutrino energies. 

motion of the bound nucleons. For 2p-2h events at low 
El mc the reconstructed energies are mostly below the 
true ones; they span the full range from nearly zero to 
i?* rue . For hihger E l J uc more and more events are re- 
conded with reconstructed energy above the true one. 
For A excitations, nearly all events would be recorded 
with reconstructed energies lower than the true ones; here 
the distribution is narrower than for the 2p-2h excita- 
tions due to the finite width of the A resonance. The 
same effect is even more pronounced for excitations of 
higher resonances, 1-pion background and DIS events 
(not shown) where the absolute cross sections, however, 
are lower. Taken altogether, for all events one observes a 
broadening with respect to the line .E* mc = El cc ; the tail 
towards lower reconstructed energies is most noticeable 
at i?* rue w 0.6 — 1 GeV. This is just the energy region 
where the MiniBooNE flux peaks; the situation will be 
similar for the T2K flux. 

The distributions for given bins of reconstructed en- 
ergy, i.e. horizontal cuts in Fig. [5j are shown in Fig. [6] as 
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FIG. 6: (Color online) 2-D density of the MiniBooNE flux- 
averaged QE-like cross section N(E tTUC , E TCC ) as function of 
true energy for fixed reconstructed energies of 0.51 ± 0.01, 
0.75 ± 0.01, 1.01 ± 0.01, 1.25 ± 0.01, and 1.51 ± 0.01 GeV. 



for understanding of experimental data. The curves are 
directly comparable with Fig. 8 in Ref. and Fig. 4 
in [21]. The difference is that our curves contain also 
the effects of pion degrees of freedom and their FSI. The 
main contribution to the observed distribution is in all 
cases given by the true-QE events, which contribute a 
prominent peak around the true energy; its broadening 
is caused by Fermi motion. In addition, the distributions 
have long tails towards larger true energies so that the 
total strength at a given true energies has sizable contri- 
butions not only from the (same) reconstructed energy, 
but also from lower ones. 

This result is in agreement with the recent analyses by 
Martini et al. [33] and Nieves et al. [23] for QE events. 
Thus, both the 2p-2h effects as well as the stuck-pion 
events lead to a shift of the reconstructed energy towards 
smaller values, or, vice versa, for a given reconstructed 
energy the true energy always lies higher than the recon- 
structed one. The effect is most pronounced at lower true 
energies and becomes smaller at higher energies. 

The cross section for zero pion events, i.e. the cross 
section for an event identified by MiniBooNE and T2K as 
QE-like scattering, as a function of reconstructed energy, 



we denote by txo,i-(-E rec ). It is experimentally obtained 
by dividing the measured event distribution (flux-folded 
cross section) by the flux at the reconstructed energy 2 . 
Each reconstructed energy contains admixtures of many 
different true energies (see Fig. [5]) so that one has 

7V(£: rcc ,£: truc )d^ ruc 

= / P(E^r\E^ m ) t t>(E^ m )a 0K (E^ B )dE^ a , (3) 

where V(E v J !C \E t J uc ) is the conditional probability den- 
sity for finding a reconstructed energy E™° for a given 
true energy i?* ruc . From Eq. ^ one can read off the 
probability density -p(^ cc |£:* ruc ): 



7>(rec|true) = P(E r v ec \El ine ) 



(4) 



1 



N(El ec ,El Iue ) 



<P{E^)a^{E^) 

The probability density 'P(recltrue) is a function of the 
reconstructed energy, while true energy is a parameter. 
It is normalized as 



P(rec|true)d££ 



1 . 



Since the neutrino flux as a function of the true en- 
ergy is fixed (i.e., an input in the GiBUU model, as 
provided by the corresponding experiments), 'P(recltrue) 
is independent of the neutrino flux and thus is one 
and the same for all experiments using identical tech- 
niques for the event identification. Fig. 11 in Ref. |17j . 
which we will not repeat here, represents the product 
7>(rec|true) aoM™)- 

The conditional probability density needed by exper- 
iment to convert the extracted event distribution, that 
depends on the reconstructed energy, into the true dis- 
tribution, that depends on the true energy, is discussed 
in detail in the appendix. 



2 Note that while in the experimental analysis the functional form 
of the flux distribution, (j>(E), is known the functional form of 
(tq w (E) is not and may indeed be different. 
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FIG. 7: (Color online) Event distribution in MiniBooNE ex- 
periment versus true (da"/d_E^ ruc , solid curves) and recon- 
structed (da I dEl cc , dash-dotted curves) neutrino energy. 



Fig. [7] highlights the influence of energy reconstruc- 
tion procedure on the QE-like events of various origin. It 
shows the QE-like event distributions (flux-folded cross 
sections) versus true (solid curves) and reconstructed 
(dash-dotted curves) energies. For true QE events both 
distributions nearly coincide. The broadening due to 
Fermi-motion, which is typical for a fixed values of the 
true energy (see Fig. 11 in Ref. [17 as well as our dis- 
cussion later in this section), has thus no visible effect. 
This in turn means, that even in the presence of Fermi 
motion, for true-QE events the reconstruction formula 
([I]) works well as long as one aims at the cross section 
measurements. This is, however, not so for the measure- 
ments of the oscillation parameters, where the neutrino 
energy has to be known on an event-by-event basis and 
Fermi-broadening thus plays a role. 

For fake QE-like events, the distributions versus true 
and reconstructed energies look very different. For 
all origins shown in Fig. [7] (2p-2h, A, higher reso- 
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nances, 1-pion background) the distributions versus re- 
constructed energy are shifted significantly to lower ener- 
gies. The largest effect comes from 2p-2h and A contri- 
bution, where the cross sections from fake QE-likc events 
amount, respectively, to 25% and 15% of the true-QE 
cross section. The overall effect for all events is that 
the distribution versus reconstructed energy is shifted by 
about 200 MeV towards lower energies with respect to 
the true energy. 

The QE-like cross section versus reconstructed energy, 
cro7i-(iv cc ), as discussed earlier, is obtained by dividing 
the measured distribution shown in Fig. [7] by the neu- 
trino flux. The flux has a peak at 0.7 GeV and falls down 
at higher energies; so for the strongly shifted event distri- 
butions (all fake QE-like) the effect of the reconstruction 
is expected to be large. 

This expectation is indeed born out. In Fig. [8] we 
show the cross sections versus true and reconstructed en- 
ergy separately for QE-like events of each origin. 3 It is 
seen that the energy-reconstruction method works well 
for the true-QE contribution, as it should. For all the 
other contributions, however, major differences arise and 
in all cases the reconstructed cross section is shifted to 
significantly lower energies. 

In one of the panels (right, 2nd from bottom) we also 
show the GiBUU cross section for all stuck-pion processes 
(Delta + higher resonances + 1-pion background + DIS). 
It should be compared with the "QE-like background" 
(open points) as obtained by MiniBooNE using the NU- 
ANCE neutrino event generator. The agreement is not 
perfect, but the shapes are very similar and absolute val- 
ues are approximately the same. While the curve ob- 
tained by GiBUU peaks at about 300 MeV, which just 
corresponds to the A excitation energy, the MiniBooNE- 
NUANCE curve peaks at about 500 MeV. The difference 
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3 Following Nieves et al. I13| we have exploited the flux uncertain- 
ties in the MiniBooNE experiment and have renormalized the 
data in Figs. [8] and [9] by 10% downwards. 
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FIG. 8: (Color online) Various contributions to the QE-like 
cross section versus real and reconstructed energies in the 
MiniBooNE experiment. The data are multiplied by a factor 
0.9. 



must be due to differences between the two generators 
used. Altogether, 'stuck-pion' events make a contribution 
to the fake QE-like cross section which is nearly as large 
as that due to 2p-2h excitations. When these events are 
subtracted from the "measured" ones, the energy shift 
due to the stuck-pion events is also approximately re- 
moved by the generator. Thus, only the energy shift 
induced by 2p-2h processes is present in the "extracted" 
MiniBooNE data. This justifies the analyses based on 
microscopic models of 2p-2h interactions (Refs. [2"31 12"^] ) 
alone (i.e., lacking the pion contributions), since these 
results are compared with the "extracted" QE data only. 

This now explains the puzzle raised at the start of this 
paper why the difference between the two data sets is so 
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large already at low energies and then decreases towards 
higher ones: this behavior is an artifact of the energy- 
reconstruction. Indeed, the reconstructed-energy curve, 
which shows the same cross section contributions vs. true 
energy, exhibits the expected behavior: it starts to be- 
come nonzero only for energies above about 0.5 GeV and 
then increases steadily with energy. 

Note that for a comparison with the data in the top 
row m Fig. |8] only the curve for the reconstructed energy 
matters. We see that this curve agrees quite well with the 
data (reduced by 10%) for energies above about 0.8 GeV, 
but overshoots them for the few experimental points at 
lower energies. We take this disagreement as evidence 
of the presence of RPA correlations. This is in line with 
results of Martini et al. [T2J (2D] and Nieves et al. [T31 
|2"T] who both find that RPA correlations lower the cross 
sections at these lower energies, but become negligible for 
energies above between 0.7—1 GeV. An inclusion of these 
correlations, which are not contained in our calculations, 
would bring our results into very good agreement with 
experiment. 

Summarizing the results of this section, we find that 
all processes other than true-QE scattering lead to a shift 
of the reconstructed energy, based on the assumption of 
true-QE, towards lower energies. For both 2p-2h and 
pion related events this shift can be quite large (several 
100 MeV). In addition the energy-dependence of the cross 
section for various reaction mechanisms plotted vs recon- 
structed energy is quite different from that plotted vs true 
energy. While the presence of 2p-2h processes explains 
the abnormally large values of the axial mass obtained 
by the MiniBooNE collaboration it is interesting to look 
for effects of deficiencies in the energy reconstruction on 
neutrino oscillation properties. This will be done in sec- 
tion [V] We note that studies along these lines have al- 
ready been undertaken in Refs. [TT1 |2"31 |2"o] . 



IV. MOMENTUM TRANSFER 
RECONSTRUCTION 

The reconstruction procedure, based on truc-QE kine- 
matics and being applied to Cherenkov QE-likc events, 
leads to distortions not only in the neutrino energy re- 
construction, but also in the Q 2 reconstruction. 

Indeed, the reconstructed Q 2 is defined as 

QL = -ml + 2E™{E tl -|^| cos 9 J , (5) 




0.5 1 1.5 0.5 1 1.5 

Q 2 (GeV 2 ) Q 2 (GeV 2 ) 



FIG. 9: (Color online) Event distribution in MiniBooNE ex- 
periment flux x dao-w/dQ 2 versus true (dashed) and recon- 
structed (solid) squared momentum transfer. The data are 
multiplied by a factor 0.9. 

As we mentioned in the Introduction, the "observed" 
Q 2 — dependence of MiniBooNE cross section can effec- 
tively be described as QE scattering with axial mass 
around 1.3 GeV. Thus, to get the MiniBooNE observed 



11 



distribution, one needs a noticeable contribution which 
falls down with Q 2 more slowly than the true-QE cross 
section obtained with a dipole form factor with Ma = 
1 GeV. In our case this is a 2p2h contribution. One 
would naively expect, that the degree of this "slowly" 
is quantified by the difference between the dipoles with 
Ma = 1.0 and 1.3 GeV. The necessity to reconstruct the 
Q 2 makes this more complicated. 

Fig. [9] shows the influence of the reconstruction proce- 
dure (JoJ) on the Q 2 distributions for the QE-like events 
of various origin. Similar to the case of neutrino en- 
ergy, for true-QE events distributions versus true and 
reconstructed energies nearly coincide. For fake events 
the reconstructed distributions (solid curves) is notice- 
ably steeper than the true ones (dashed curves). This 
leads to the same effect for all QE-like events. Thus, the 
reconstruction procedure |5| makes the Q 2 -distribution 
look steeper, which in turn means, that the distribution 
of 2p-2h contribution versus true Q 2 should be event flat- 
ter than the naive expectation. 

Within the 2p-2h model employed in this paper, for 
"2p2h+QE" curve (right top panel) the agreement of 
the reconstructed curve with the "MiniBooNE extracted" 
data is not perfect. For lower Q 2 the calculated curves are 
higher than the data; this is the region where RPA effects 
should bring them down [HHT]. For Q 2 > 0.35 GeV 
our reconstructed curve is steeper than the data. For 
all events ("MiniBooNE measured", top left panel) the 
differences are larger, which is due to the different treat- 
ment of stuck-pion events in the GiBUU and NUANCE 
generators. 



both studies the QE-like events were used and neutrino 
energy was reconstructed assuming QE kinematics, ac- 
cording to Eq. ([TJ. Thus, the effects of energy reconstruc- 
tion must be similar to those in MiniBooNE experiment. 
In this section we show, how this influence the extraction 
of oscillation parameters. 
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FIG. 10: (Color online) QE-like event distribution for the 
original (thin curves) and oscillated (thick curves) QE-like 
event distribution in T2K flux for muon neutrino survival 
measurements. The same events are shown versus true 
(dashed) and reconstructed (solid) neutrino energy. 
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V. EFFECTS ON DETERMINATIONS OF THE 
OSCILLATION PROBABILITY 

T2K experiment has recently reported the first ex- 
perimental observation of electron-neutrino appearance 
from a muon neutrino beam |27j and investigated muon- 
neutrino disappearance |10] with an off-axis beam. In 



FIG. 11: (Color online) The same as Fig. 10 but for elec- 
tron neutrino appearance measurements. The oscillated event 
curves have been multiplied by a factor of 10 to enhance the 
visibility of their difference. 



Fig. 10 shows the actual oscillation signal plotted as 
a function of both the true (dashed curves) and the re- 
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constructed (solid curves) energy for the muon-neutrino 
survival events. Fig. [TT] show the oscillation signal for 
electron-neutrino appearance events under assumption 
<5cp = 0. Thin curves show the QE-like event distri- 
bution for the original T2K flux [10], while thick curves 
that for the oscillated flux calculated with the oscillation 
parameters 



Am 2 3 = 2.5 • 1(T 3 eV 2 



and sin2#oc, = 1.0 



for the muon-neutrino survival signal, and 

Am 2 ,! = 7.6 • 10~ 5 eV 2 , 
sin 26»i2 = 0.927, sin 29 13 = 0.316 



(6) 



(7) 



for the electron-neutrino appearance signal. 

In the disappearance signal the oscillation minimum 
is significantly affected: plotted as a function of recon- 
structed energy the minimum is smeared out and shifted 
to a higher energy (by about 50 MeV). Also the first 
maximum is similarly affected. Correspondingly, in the 
electron-neutrino appearance experiment the maximum 
vs the true energy lies higher by about 20%. 

To make these points even more visible and to demon- 
strate their effect on the oscillation probability we show 



in Figs. 12 and 13 also the oscillation probability defined 
as the ratio of the oscillated to the initial event distribu- 
tions. Again the minimum in the curve vs reconstructed 
energy is smeared out and and shifted by about 50 MeV 
towards higher energies. The same is true for the first 
maximum in the electron-neutrino appearance probabil- 
ity. 

The effect here is, however, not so large as for the event 
distribution. This is because the energies for the original 



flux (thin solid curves in Figs. 10 and 11) are shifted in 
a similar way as those for the oscillated flux (thick solid 
curves) . 

It is interesting to compare the uncertainties in the 
oscillation signal that are caused by the energy recon- 
struction with the sensitivity of the oscillation to a CP- 



violating phase. This latter is shown in Fig. 14 It is 
clearly seen that the sensitivity to a CP violating phase 
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FIG. 12: (Color online) Muon neutrino survival probability 
versus true and reconstructed energy for the T2K flux. 
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FIG. 13: (Color online) The same as Fig. 12 but for electron 
neutrino appearance probability. 

is comparable to the uncertainty caused by the energy 
reconstruction which can be seen in Fig. |11| 

VI. SUMMARY AND CONCLUSIONS 

The neutrino energy reconstruction using kinematics of 
QE scattering depends on a clean identification of true- 
QE events. In nuclear targets this, however, becomes 
difficult. Both new primary reaction mechanisms and 
final state interactions may cause a misidentification. 

A recent example for the errors introduced by such 
a misidentification is provided by the energy-separated 
cross sections for QE-like and QE scattering obtained by 
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FIG. 14: (Color online) QE-like event distribution for the os- 
cillated QE-like event distribution in the T2K flux for electron 
neutrino appearance measurements for various CP violation 
phases. The oscillated event curves have been multiplied by 
a factor of 10 to enhance the visibility of their differences. 

the MiniBooNE experiment. Here the extracted QE scat- 
tering distributions could only be fitted by using an axial 
mass of Ma = 1.3 GeV and this result has lead to a dis- 
cussion of possible in-medium effects and to a multitude 
of partly contradictory explanations. Now the general 
consensus is that this result was obtained because the 
true QE scattering mechanism was not correctly identi- 
fied and, therefore, the extraction of Ma was made using 
the wrong theoretical framework. As a consequence of 
this incorrect identification not only the extracted Ma 
value was incorrect, but also the QE cross sections as a 
function of reconstructed energy were found to be wrong. 

So far, only the QE scattering cross sections as ex- 
tracted by MiniBooNE have been analyzed [23 El . This 
extraction involved some model dependence since the so- 
called stuck-pion events were removed from the experi- 
mental QE-like data set by using a tuned version of the 
generator NUANCE. Motivated by this fact we have, in 
this paper, analyzed also the full QE-like cross sections 
in one consistent model, GiBUU, which contained all the 
pion production events, true QE scattering and the 2p-2h 
excitations. The main result is that all other processes, 
not just 2p-2h, lead to a significant downward shift of the 



reconstructed energy. For the MiniBooNE flux this shift 
amounts to about 500 MeV for the 2p-2h excitations (at 
the peak of their event distribution), to about 500-600 
MeV for pionic processes involving the A resonance and 
pion background contributions and to a shift of about 1 
GeV for higher resonances and DIS. On the other hand, 
for true QE scattering the energy reconstruction works 
quite well (by construction). For the full QE-like event 
that involves all these reaction mechanisms the shift in 
energies then amounts to about 200 MeV. 

Even larger is the effect on the energy-separated cross 
sections where the QE-like and the extracted QE cross 
sections have a dependence on energy that is very differ- 
ent for the true and the reconstructed energy. We have 
shown in this paper that the analysis involving all reac- 
tion mechanisms can explain both the extracted QE and 
the measured QE-like energy separated cross section. It 
also explains naturally why the difference between QE- 
like and QE-extracted events is quite large already at low 
energies and then becomes smaller toward larger neutrino 
energies, contrary to expectation. This unexpected be- 
havior is again due to errors in the energy reconstruction. 

We have, furthermore, shown in this paper how these 
uncertainties affect the analysis and the extraction of os- 
cillation parameters from the ongoing T2K experiment. 
Most important here is the result that the uncertainties 
necessarily connected with the energy reconstruction are 
as large as the expected sensitivity of the oscillation re- 
sult on the CP invariance violating phase. 

Finally, in an appendix, we have given the transforma- 
tion from reconstructed to true energies in a form that 
may be useful for experimental analyses. 
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Appendix A: Energy reconstruction 

In an experiment only the flux averaged cross section 
is directly measured and by the reconstruction method 
the product 4>(E Tec )a 0l7 (E Tec ) can be obtained. The chal- 
lenge then is to obtain the true event rate from this anal- 
ysis. 

The true distribution can be obtained by integrating 
out the reconstructed energies in Fig. [3] 



^{ET C )^{ET C ) ■ (Al) 



This true event distribution can be related to the recon- 
structed one by 

0(^ ruc )a o .(i?* rue ) = 



where V (E^^E^) is the probability density to find a 
true energy El IUC in a distribution of events all having 
the same reconstructed energy El cc . 

From this equation P(true|rec) can be read off. It is 
given by 



-P(true|rec) = T(E" uc \El cc ) 



(A3) 



1 AT f 77TCC 77'trUC\ 

Here <f)(E vcc ) stands for the true flux, read off at the 
reconstructed energy, and the QE-like cross section ver- 
sus reconstructed energy, a 7r(El cc ), is defined by Eq. ([3| 
and shown in Fig. |J In Eq. ((A3]) ££ oc is a parameter 
and £?* rue a variable, so that the probability density is 
normalized as 



7 3 (true|rec)d^ n 



1. 



The distribution 'P(truelrec) is of primary interest for 
neutrino experiments. For fixed values of true and recon- 
structed energy, its connection with the 'P(recltrue) can 
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FIG. 15: (Color online) Conditional probability densities 
"P(true|rec) of finding a true energy _E* ruc for fixed recon- 
structed energies of 0.25 ± 0.01, 0.35 ± 0.01, 0.65 ± 0.01, 
0.95 ± 0.01, 1.05 ± 0.01, 1.25 ± 0.01, and 1.35 ± 0.01 GeV. 



(A2) The fits are according to Eq. (A5l 



TABLE I: Fit of the :P(true|rec) according to Eq. |A5 
) 
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be obtained by comparing Eqs. ^ and (A3), 

c/>(£r c K.(i?r c ) 



"P(true|rec) 
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V (rec| true) (A4) 
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The probability densities (A3) are shown in Fig. 15 for 



various reconstructed energies. We have found that in 
the energy region 0.35 GeV < E™ c < 1.25 GeV they 
can be described quite well by a sum of two lognormal 
distributions: 



P(true|rec) 



1 



V2tt 

1 — Hi 

xa 2 



X<J\ 



exp 



exp 



(ln(x) — ln(mi)) 2 

(ln(x) — ln(m2)) 2x 1 
2a! 



(A5) 

Here x = i5* ruc is a variable, m\ = E™ c is the center of 
the first peak, which (as discussed above) mainly comes 
form Fermi motion in truc-QE scattering. The second 
much broader peak is due to the 2p-2h and stuck pion 



QE-like processes. At low E r v ec it is centered at notice- 
ably higher energies then the first one, with increasing 
E r v ec it shifts closer and closer to the first peak and even- 
tually to the left of the first peak. Parameters n±, ci, 
m.2, and a 2 are fitted. In Fig. [l~5"| the best fit is shown by 



the solid line, the first term in Eq. (A5) by the dashed 



line and the second one by the dash-dotted line. 

For other energies the fit parameters are given in Ta- 
bic u 

For lower, E™ c < 0.25 GeV, and higher reconstructed 
energies, E'J 10 > 1.35 GeV, the shape of the probability 



distribution differs from the functional form (A5). This 
is also illustrated in Fig. [15) 
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